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A \eo-demenaonal mantel 10 devcribe clgt-eMects of planet cuculer dec electrodes has bere devetuped T 
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erturgonad colle ston 0 ured te ducienss the 
chrsusmperoecins capenenadd era dared. 


wo deme vonad 
equemens. The equetune tor the ueadeuon of a 
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On pes es punet une methade bedamenusanetie pour decre ies effets d' artic an arvese des Clectredes 6 desques 


cuculewes 


plant (ino widoe une dnpounen arihogunalt pour let cquatone deflerenuciies parueiies correspemdantes. On o denve irs 
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fntredurtion 

One common simphGcation ia the theoretical 
treatment of electrode processes of planar desc 
electrodes 13 the assumptun of renv-iehnite hnear 
diffusion. Thes results in the reduction of the 
three-chmensional sysiem to | one-dimensional 
mudel, where the distance + from the electrode 
surface is the only courdnate to be comsdered. 
Consequently, the Uwee-demensonal form of Fick's 


second law : 
éc ae oe ae 
“ a poe 
reduces to the sumpler differential equation 
ae ae 
Beas 
and the mathematical end computational efforts 
decrease . The essumptun of hacer 
deffusson t0 2 planar electrode angi unply an i- 


Rane dic without edges. For real systems, how- 
ever, thes of course ts not possible. 

Thus, measurements of fiente diac electrodes 
should devuste from the simple one-demenvonal 
theory. In fact, t has been shown for delflerent 
eleciroanalytical techreques such as chronuem- 
perometry and chronopatentiometry (1) thet the 
dofluswn current 19 bogher than expected from the 
theoreucal calculations based on eq. [2], 1.¢. the 
Cottrel equation (2) fur chronoumperometry. For 
cychy voltammugrams (3) the peak separation 1 
tasger than the theoretical values geven by Nochol- 
son and Sham (4). A number of papers have 
accounted. for these deviations as being duc to 
addnvonal deffeuon from dwechons other then 
perpendicular to the dinc electrode (ace ¢.g. refs. 5, 
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6}. Ina system with several coupled charge transfer 
and chermcal reactions dramatic changes of the 
overall voltammogram have been attributed to 
edge effects (7). 

Several attempts have been made to calculate the 
effects of edge difluvan and to simulate model 
SySIEMs without assuming onedhmensional dilTu- 
swon (5, 6). Recently, Hemze published a two 
dimensional model bared on the firute ddference 
methud (8, 9), while Oldham (10) and Aoks and 
Osteryoung (13) developed analytical soluuons fur 
the umple reversible mechanism. 

These calculatiuns showed, in agreement with 
experiments, that edge dvfluvion rates mcrease with 
sncreasing ume and decreayng radas of the disc 
electrode The deviation frum pure Cottrell behav- 
jor 13 claemed to be between 9 and 10% for typecal 
conditions. diffusson coefhcrent = 10 * cm/s, clec- 
trode radars 0.1 cm, electrolyss ome 2.655 (5, 9-11). 
All methods presented thus far, however, demon- 
Strate results for unly the somplest mechanisms 
(e g. ceverseble electron transfer case). Extensions 
apparently requere extensive compulatvonal cifurt. 

Orthogonal collocation has been shown tu pro. 
vale a powerful method to simulate responses of 
several clecuroanalytical techniques to a varicty of 
mechaersms (12-18) Up to nuw all these calcula. 
trons have been performed using one dimensumnal 
deffuscon models. We want tu repurt the consiruc- 
thon of a two-dimensional model for the diffusion to 
and from a planar dive electrode and the dixcretiza- 

ton of the resulteng second order partial differential 
equations tu give a set of semultancous first order 
ordinary ddlerential equations. 

The companson of the results of this new model 
with those already presented is the subsect of © 


esoeenl do recherches do Cosas 


4 


TSS eee - 


SPEER AnD pons 


Separate repurt now in preparation. Extensions to 
other mechanisms and electroanalyical expen- 
ments will also be consuereu. 
The two-dimrnsinagl ensdet 

Consider a planar, circular disc of conducting 
Material ambeukled into an mvulator (Fig. 1). The 
radius of the disc is re. The radius 7, of the 
insulator may go to infinity. From the symmetry of 
the problem, follows that we may use two 
Courdinues £ (perpendicular to the disc) aad 1 
(along the sadws vector of the disc). We ase 
fequued to calculate the concentration ¢ of all 
Species as a function of 1,7 and the tune ft, where c 
és given by second order parual ial equa- 
tions of the type (9 


oc ae bade. de 
Ox [35 7+ F5]-20 
and fic) is determined by the kinetics of the system 
investigated. J 1s the diffusion coefficient. 
Let us consider an EC-mechanism (4). There the 
Gifferenual equations descriting the 
of species A, B, and C are gsven by [5}-{7]. 


4 Amie atic 

(sh tem Dy ste? Say Be 

a md ere ed ee 
i" ra + hcp 


For & = 0 we deal with a simple reversible electron 
transfer. The electrode reaction m eq. 14) is drawn 
8s an ordain bul the eatensiun to reductions is 
fnmedtate. Fur simplification we may assume thet 
the deffusion coeffictents are equal: 


(} HDD, 


Derumetnc expenment, where species A, the sub- 
Sirate, ts present at a concemtration 
2 at the beginning of the expenment. The product 


4359 


Pro. 1. Geomaey of the twe-deeseninnsl ented. 


Species B and C may not be present af ¢ = @ (initial 
Condhtions (9), | 104) < 


$=6,120,°26 
I} clr. Mack 
(10) chs, 7, Mocs, 4, HOO 


t>0, 5-00, 720 

(U} egies, thee 

U2) gh, 7, he cee, 7, he ® 

The same condition may hold for a distance 
from the electrode radually: 7 = 
>6,226, 7-0 

11) ele, 9. eek 

[ta] cule. @, oes, ©, 8 

At the electrode surface the Nernet equation has te 
be fulfilled for c, and re (eq. [14a). M we step toe 
Putentual E ot least 200m V past the formal potenual 
EX~ the concentrauon of the substrate A at the 


és wamedstely forced to zero and [14a] 
May be replaced ty | 146). Additional condstions at 


0384 cam. §. comm 
the electrode surface are given by the Bux eque- 
tions (15) and (16). 
(> O, ce @ OSrS% 

(0.0 aF ce 
{14a} AO wane [FU | 
(146] (0,7, 8 

dea 


ee 
ua (%)_--(F). 
dee 
—_— e 
mt (5° 
Ata = 6, but outside the efectroactive area(s > re). 
no Species may undergo electron transfer and we 


have boundary condwion {17}. 
t>@, 1520, Sree 

dea) (22) 2 (2) 
mm =). (S)_. G)_. ; 
Because of the symmetry of the problem there may 
be no act fun of any species across the x-a21s at the 


comer of the electrode (7 = 0. 
s>0,220,7=6 


we (S)-(F-G)- 
Equations (4}-{18} describe the two-dimensional 
enadel system for an EC-mechanisat ender chrono 


2 = O inside FS %e) and outsde (7 > 1) the 
electrode, it is convenient to describe the whule 
system by two coaceniration 

which are defined for 05 7 S ro andr < 7 
respectively. To avoid discomimsities at 7 = fo WE 
have to apply some additonal “boundary” condi- 
trons for all specees: 


vor. as, 19 
In order to meet these requirements we use the 
following iransformaiion ¢quations: 


(2) Xealh 
(22) c* mek 
(23a) R= rir, 
(239} R= - 1 KM - 10) 
(24) eeh-ryD 
(25) Te Debi 
where X is the dimensionless distance from the 
electrode, L is the distance from the electrode ia 
the x direction defined so that at X = L no changes 
ia concentration take place during the time of the 
expenment, c* is the dimenwonless concentration, 
R and R’ are the dimensionless distances ia the 7 
director inside and outsule the electrode disc 
respectively, Mf is the distance im the 7 dwecuon 
defined oa the same manner as L, and T° ts the 
dimemsioniess simulauion tne. 

Application of eqs. 1211-125] and the subdivision 
into inner snd outer functions c, and c, of eqs. 
{4}-(18) yield the dimensionless mode! [26}-133} 
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or’ ox? oR aR? 
och a'eg, . 1 dek ack, 
ON Spee ae te on tone 
och Bed, 1 ded, , Teh, 
Om pel ag te oe tone teh 
where 
(24 P= ruL? 
and 


def, d'chs w = decks 
DH SF oo teaver 


aad, eee ee . 
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rane ae T'>@, X00, 0c A's 1 
a cI M - ry dels éc8, ae 
with the initial condtteons iad ox ). (3 )° te 


T'=0, OS X51, OS R<1, Oc At 

1394) BUX RD LAX RO I 

(35) RX, R= RIX, RM — ChAT, RO 
= c8AX. 8.8 

and boundary conditions 

T'>0, Y=l,O<As1 

O64 chILA TIO! 

D7 ch CL A. TI — 801, &, TI e 


T>0, X=, 0sK's1 
38} cht are! 
OF ch RT| cha Pee 
T>6, 0s Xs1, R'=1 
(Of 2AX, 1, TOI 
(40) BAX, LTV CALI OO 


T>@, X=0,0<a<1 


CLO, RT) [ak 
Wal 0. ere (are -| 


(2d) 20,8, TI = 6 


or (ar).7- (ar). 


ack, 
wa (FF) 78 


T>@ es Xs), R08 

m (SL (SEL (GP -* 
T>@, Os X55, Ro1(Ke® 

I RU = 8 hd. 8,19 


HO CRUX, T= chAX.8, TD 
Me A) 


om (SL (SEL 
on (StL vee 
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tions of species A, B, ond C respectively 
function of X, AUR, and T’. = 


153) che, #, 8) Ay pA Ag eH... Hag gg ge OP te 


Hees ttrOOO! + Dan opt OOOt Ole sag ge brtysort 


eae paula iad +p ech + DagtOr! + ggt'e! +... Hd ons tet ty? 
‘The second derivative of ¢ with respect tor in only ont + He eae 
dohiel i i Ne ere ee or a eb db oh, 7 
om aren farther —— = § — 
cha ea as PO or Par revi OF ae 
ae) (Se +p 2h _ me 
“i (F), (=), Pont ere (s4) ctsesne SS nasties 
wines been pointed ext (17) that orthogonal (221 oh ag eh, Fe ei 
: om 8 ox! oR’ +Vp OR i.e. we mi the i 
collocation may be applied ealy to dimeasionicss MF art, ina oad’? whose aaedspeaasn ccahcunte 
syviems which have been converted from som +0 5 tech the 2,,. The polynomial is of order N, + 2 and N, 
infinite models to models in the range @ 00 1. oR’ + 2 with respect to # andr, respectively. The 
nero : aA, Sa peg Nees So 
j 
~N _ 
ae 
. 
- et = a wee ee 
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points in the s and r directions, respectively 
Rea med 

IH) cine FE maeet 

The concentrations of the collocation points may 

be the elements of a concemranon matrix I, 

(Tg Mer ABE 

which is calculated by mubkiplication of a matrix 

X, 5. contasmng the 2/°' as elemeals, a mais A,,, 

Comtameag the elements A,,, and & mainn RZ,, 

Contamung the elements ri '. The superscnpt T 

denotes a transposed mains. 

Now we can calculate the vst and second 
denwvatives of the concentration af the collocation 
ports with respect to s and ¢ respectively: 

(SM Peg Xfy Ags RE, 
$8) FAB RSA Rls 
(9 Me By As BES 
Te hy A 
where the saperscnpt s(r)} denotes the firs! and 
accordengly s0(r7) characterises the second derive- 


bove wh respect to #17). 
From 0. (56) 0 follows thes 


I Ay OX) Fe (RL! 

After iqnerting (61) into [57] we arvive at 

et Cty RO Kes Oya A Te 

where the fertt deviveteve of the concentretion with 
rerpect to sm enprened as a functinn of the 
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From (61] and (99 we get 

Pe MAT MP ar eee 
Muy AML, 8,37" 
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Analogous from (61) and (60) 

fee) Fe (BT LO Ny Be 

We mav cakulue the masices containing the 

dervvatives of the cancentrauen af the collocation 

prents from the conceatrawen mains ence the 


muirices A,,, B,, Cs, and R, have been deter- 
auned. Rewsung to wammetion express we 
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{e7) x) = . ps Ascly e. 0) 
ay 
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(es) =), - bs Bisa fast) 
ae mer 
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Pee) 
and 
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ae 
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As has been pointed out by Whiting and Carr (12) 
far the one-kumenssonal case. the matrices A,,, 
By. Ci». and DO, depend solely on the type of 
polynumal used as tral function and may be 
cakulated from the zerues of this pulynonual. This 
has to be done caly once. Using the values of the 
elements of those malswes we are able to express 
the defferentis! quatients of c with respect tos andr 
at ail collocatron punts. 

For bnear ddfuson problems it has been shown 
thal Legemire polynomials are most susable (12). 
Pons (18) used Cheby.heff pulynommals for diffe- 
son wader sphencal geumciry. It may be nuted, 
however, that the chusce of the trial function dues 
won have a large effect on the accuracy of the 
soluteuns for Ingher urdet approxsmations. la gen- 
ge ee 
symmetry of the problem. 

ta thes work, Legendre polynomials have been 
weed far deacretizauun along the 1-courdnste. The 
collocation in the drectwmn of the ¢-coordhanie has 
been performed using diferent types of functuns 
wnchodeng Legendre, Cheby vheff, aad pulynunuals 


We can subtinute the tummauon expressions 
(672~-(701 for che parted quuhents m eqs. 
(20]}~{ 28) and {30}-(32} Thes yrehds expresaouns fur 
the tome dervative of the concentration of species 
A, B, end C at alt the collocation pots 23 8 
function of the concentrations at the cullucaton 
ports. we subsisute the boundary conduwes 
[301-{46], (4201 (i.e. we assume a step to a potenual 
Gas behind E*) and |47}-(49), we arrive at 


awh, 
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r feractur 0.774 Exssatht%e T+ Fnac. &. n} ~ect, 
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+ [eserix.0 794 En 2.268% 1TH ¥ eactita m. n} rock 

ra hk. 6 la sct00. R714 Raat “Smet Re n| 

” [riseta 079+ Raat S Fisch Ri | 

ns) wh =P lmsetate. R794 "S mac lts Re, nj 
+ [riachdt. 0794S Fisctlte R ry) - et 

TEU 0 (aactae nce S acts, n| 


+ (risctase, 714 risctdte Ri | tects 
where 
TM Byn Set Ree 
end 
(781 Fi,* 


Cie + Dy, 
VO OR * 
‘The prmed terms, except T° nd 8’, correspond to the outer functions. j é 

* i) 
+ ry : Ie vaknowns, 1.€. 1%, 0, TD. 20. TD ChK Trete, Theme’ Th 
CHA 1. TD, CRAM, 0. TC BAN OT, CHAN, ©, T)bor alt Xi; CHAO, Re, Tend 8,0, Ry, T°) for all 
boundary 


Ry. cL HO, mrs ch4@, mT), RAO, Ry, T°) for all Rj which may be computed from the 


Comiason (43) yseids 
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™ ¥ Achy Me, T= — Sac M1) 


ciectde teens CE ne 190% (571, ent 420) end safving for the concentration of wpacias Bot the 
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Condition (44] is changed to 
"S auchtt Ra M20 
aia ope, aot tnstinsion a8 7% 
(SI) cBMOL A, Den ge Faget RTD 
Sa anadagy it may be derived from [45] and (40) or [41] 
(1 MATIZ [Aenea S auetan 7] 


-_— 


Roveaéion 7 
REIS) G>:, 
DTIc r°3 
Unagmounso¢ 
Sustificetto: _ 
ee 
Distrtout tons 


| i a 
cht. Ri. Te -2¥ Ay thX, BT) 


EEE 


1 aon 2 
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‘The condition [46] gives, for species A, B, and C 
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which may be expanded and transformed using [40] and (47) 
199) Lemer.nyer— VP Chak e8 te, 1, 174 Conger ®, TY 
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Conditions [46] and (52) (146] and (S2p ere weated simiarty using the corresponding boundary ex- 
pressiuns. Thes, 


oe | 
199) chIM LT hd. T Ie z (¥ Metts. Ras 71+ F CaVP Cite an | 
and 
r oa e 
I) CAML TIO ctdtne.TI= 2 (¥ Meh Re nS Cua VP Cuctlty RETY 
Finally, the concentrations of species A af R @ Gand all X, nay be calculsted by introducing (90) into [85}: 
195) chX.8, re-z [Sacnu.a.7 
> Cason VO Cis tA, Ris TY 4 Casas WP Gipos } 
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PY We. 2 [s Gentte. &. ne¥ Cuens VP Cette Re n| 
and that 
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Now we any substions the unknown terme in ome. [71}-{76). Prom (71), using (90) and (95) it fellows 
ets 
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Combination of [72] with (80), (93), and (96) gives 
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The outer functions for all three species are deri 
CT in the outer domaia is found from [74], [82]. 
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the same scheme. The time derivative of 


derived following 
‘and (90). and ia analogy the correspondeng equations far 


species B or C may be calculnted from (75}, (83). and (93) or [76), [84], and (9). 
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Now the set of six second order partial differential 
equations (26}-(28] and (30}-{32) has been trans- 
formed imo a set of 3N,(N, + Ny) ssmultancous 
ortenary differential equations |98), re (103), 
(104), Nar ep al Such systems may be casdy 


2d C af the colacatiue puts (X,,R,) and(X,, Ry). 
We describe the concemtratwon profiles winch de- 
velop in from of the a by the concentra- 
tions at the collocation pow: 

We note the sacar say fhe ferent 
dexcribing the tims derivatives of the Concenwa- 
tonne. This ie one of the advantages of the orthog- 


i 


onal collocation meihod, for if the theory hes been 
derived Gace it is very easy to express & im 
computer laaguage. Furthermore, # more compli- 
cated mechanisms are to be semulated, only arnur 
changes in the formulse have to be performed. 
Also, the extension of the theory to other experi- 
mental techraques, i.e. other boundary coaditivas, 
és casy. This will be shown im another paper of this 
serves. 

We sho note that the i of the 

expressions descriteng species C is aut essential to 
calculate the concentration profiles of A sad 8, lor 
C te not electreactive aad dues not react lo or with 


an clactroactive mtermedine (16). However, fone 


eae ent, 
SSS ee 
tote 


desires, for example, to calculate the spectraciec- 
trochemical texpunse (14) of species C one is 
interested in the concentration pruties of C us well. 


where A és the electroactive of the electrode. 
‘This umphes tha the ft of the bette wert 


This is cleasty not true im the case of a two 


dimensional mudel consalering effects. 
Stead, the current is given by ~~ bed 


110% i= nFD [lr (=). « 


where the r-dependeat flux through an infnitesi 
mally small part 2arde of the electrode i ianegrated 
we the electroactive section 0-1, of the >- 
Using the transformation equations [21}-123) 
(251. and the definuion of p (eq. Prag dy pig 


ante a a [ia (SB) 


= mract 2 VF fin (SE) a 
‘The chroncamperometric “constant” {Vj ie 


given by 
WU tie race VBVET 


fin (ar). 


Abernatively, we press 
terms ofthe Cottrel comsient™ ¢2) Ot 


(012) vi - FAckVD 
te Vee 


ad 


Since the 

xi Socotinn poe ot CBP) ae given 
irtegral exactly by a quadrature fone 
Pian ta forma ofthe type 
ttf tomine= San 
. e 

Warze the Q) are elements of « vector depending on 
Thus, from (67) 


"18%, 
om [et 
ot 

~ 30 [a Seucton.a.r] 
Expanding the sum over j and 


bouadary Conditions (36 nad (20), [15] cRenges 


we f R SEs an a 
* [Petra +h x Aetihy &, n| 


The Geussion 
are equal to ero, ee 198 eae Ne #2 


wn fia (HE) z a 


" [Marner m SF aueta Rr, n| 


which is 8 conventions form for 
imegral, computation of the 
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Intredactiea 
The first dinuclear chodium complex of the 
“A-frame” (1) [Rody-S> 


geometry devignaed 
(CU) Adpm),) (dpe is Ph, PCH, PPh, and dam is 
Ph, ASCH, AsPh,) was descnded in 1977 (1). The 
synthesis aad siructure of complexes of 
thodnim of indium has since received muck attea- 
tion (2-16). For complexes of rhodsent|), several 
instances of co-ordination of two-electron donur 
racecar CO, SO,) at the site dridging the 
metal centers (2) have been reported (2. 3, 6, 9, 
#210). The complex cation, Laretrpeg 
(épm,!*, veobagiel 2 1:1 kabsle, symmetrical (, * 
am) complen with diphenylacetylene (6),? sad 
semalar neutral complexes have also recently been 
described (11). Co-ordination of an alkyne at the 
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